arXiv:hep-ph/9507407vl 25 Jul 1995 


hep-ph/95U7l07 

RAL-95-035 

OUTP-95-13P 


DO 'll) (4040), 'll) (4160) SIGNAL HYBRID CHARMONIUM?f] 


F E Close 

Rutherford Appleton Laboratory 
Chilton, Didcot, Oxon 0X11 OQX, U.K. 
fec@v2.rl.ac.uk 
and 

P R Page 

Theoretical Physics 

1 Keble Road, Oxford, 0X1 3NP, U.K. 
p.page@ph.ox.ac.uk 


Abstract 

We suggest that (4040) and U (4160) are strong mixtures of ground state hybrid 
charmonium at ~ 4.1 GeV and the 'f{3S) of conventional charmonium. The T® ® , 
masses and total widths of the ^'(4:040) and ^’(4160) are in accord with this hypothesis. 
Their hadronic decays are predicted to be dominated by the iplSS) component and 
hence are correlated. In particular we find a spin counting relation r(4160 —> DgD*) ~ 
4r(4040 —> DgDs) due to their common '^{38) component. For D and D* production, 
using ^/)(4040) branching ratios as input, we predict that the decay pattern of the 
■0(4160) will be very different from that of the 0(4040). These predictions may be 
tested in historical data from SPEAR, BES or at future Tau-Charm Factories. 


^Work supported in part by the European Community Human Mobility program “Eurodafne” contract 
CHRX-CT92-0026 and the University of Cape Town. 




The QCD sector of the Standard Model cannot be regarded as hnally established while 
the qnestion of the existence of glneballs and hybrids (hadrons where the glnonic modes 
are excited in the presence of qnarks) remains nnresolved. This is now becoming rather 
critical. On the theory side a consensns on spectroscopy and dynamics is emerging from 
lattice QCD and model simulations thereof whereas on the experimental side, the increasing 
concentration of resources at the high energy frontier threatens to leave a hole in this 
essential area of the Standard Model. 

We have initiated a programme to evaluate the best opportunities in this area and have 
identihed possible signals worthy of further study 0, H, 0. In this letter we note that 
e’''e“ annihilation in the vicinity of charm threshold, 4-4.5 GeV in offers special 

opportunities: the lowest hybrid charmonia “ifc” are predicted to exist at 4.1 - 4.2 GeV 
IQ, Q], just above charm threshold where the conventional cc are rather well understood 
and where there is a well known experimental anomaly which may be due to excitation of 
these hybrid states Q. As a result of our recent investigations into the dynamics of the 
quark - flux-tube system [Q, Q] we present new arguments that support this suggestion and 
propose further experimental tests. In particular we shall argue that ^’(4040) and 'ip{4160) 
are strong mixtures of '0(3*5') and hybrid charmonium. 

region of the following those of ref. [Q] show 

One of us has shown |Q| that, within the conventional charmonium picture, the hadronic 
width and branching ratios of "0(4040) and the narrowness of the -0(4415) constrain the 
-0(4040) to be the 35 (as opposed to 2D) and the 0(4415) to be 4S (or possibly 5S) states; 
this independently supports the spectroscopic assignments of many potential models. The 
0(4160) hadronic width is consistent with it being any of 35,45 or 2D Q] and hence, 
by elimination, one is tempted to deduce the 2D or possibly the 45 Q] assignment for 
it. However, r®®(4160) ~ r®®(4040) which is inconsistent with any of these pictures. The 
simplest explanation is that these states are roughly 50 : 50 mixtures of 0(35) and an 
“inert” state that is essentially decoupled from as originally suggested by Ono 0. 

The unresolved problem is what “inert” state the 0(35) is mixing with: dynamics 
within standard spectroscopy appear unable to generate such a mixing consistently. Mixing 
between 0(35) —0(2D) via a tensor force appears to be inadequate given the known limited 
mixing for the 25 — ID states 0(3685; 3770) Q, ||, 0. Strong mixing between cc and DD 
coupled channels is incompatible with the observed narrow widths (of order tens of MeV) 
of the 0(4040), 0(4160) 0. The maximum theoretical mixing of 0(35) — 0(2D) generated 
in the literature, to the best 

of our knowledge, is some 10% Q, Q] in amplitude. 

However, a mixing of 50% can arise naturally if there is mass degeneracy between two 
“primitive” states. The 0(35) and 0(2D) are not expected to be degenerate enough for 
large mixing .0 However, a new candidate for an “inert” state degenerate with 0(35) has 
recently emerged. 

Numerical solutions of flux-tube charmonium spectroscopy predict that the lightest 
hybrid charmonium states are in the 4.1 - 4.2 GeV region Q, while QGD inspired potential 
models with long 

range linear behaviour uniformly predict that 0(35) is also in the narrow range 4.10 to 

^ A coulomb + linear potential gives ||^ a 2D-3S mass difference of 79 - 87 MeV for string tension 
b = 0.18 ± 0.02 GeV'^, a charm quark mass iric = 1.5 ± 0.3 GeV and as = 0.4. The ranges of b and TOc 
are the maximal ones consistent with spectroscopy. The only way of reducing the 2D-3S mass difference 
signihcantly appears to be by reducing as; nonetheless it is still ~ 30 MeV even in the extreme limit 
as 0. 
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4.12 GeV p, 1^, 1^. This near degeneracy is also snpported by lattice stndies of cc which 
predict that '0(35') and are within 30 MeV of each other [Q and, again, in the 4.1 - 
4.2 GeV region. Gonseqnently it is probable that, within their widths, there will be mass 
degeneracy of 0(35') and He, leading to strong mixings and splitting of the eigenvalnes. If 
snch a degeneracy ocenrs one immediately expects that the physical eigenstates will tend 
to be 

V>± - ± He) (1) 

which we shall identify as 0_ = 0(4040) and 0+ = '0(4160). 

The He component in eq.(l) will be “inert” for the following reason. A characteristic 
feature of hybrid mesons is the prediction that their decays to ground state mesons are 
suppressed and that the dominant coupling is to excited states [^, |^, Q, in particular 
DD** for which the threshold is ~ 4.3 GeV. Hence this pathway will be closed for ^(dOdO) 
and '0(4160) and consequently their hadronic decays will be driven by the 0(35') component. 
Thus the He component immediately satisfies the “inert” criterion for the additional piece 
and is partly responsible for the relatively narrow widths of these states. 

Both the dominant production in e+e“ annihilation and the prominent hadron decays 
will be driven by the 0(35') component, leading to intimate relationships between the 
properties of the two eigenstates. First, this naturally explains their leptonic widths 

which is consistent with the fact that potential models uniformly predict a value for r®®(35') 
that is essentially a factor of two larger than the data in both cases. 

We now consider the implications for hadronic decays with particular reference to our 
recent dynamical analyses B 1,1- 

The remarkable branching ratios (with phase space removed) 

0(4040) ^ DD ; DD* ; D*D* ~ 1 : 20 : 640 (2) 


differ considerably from those expected for a simple S-wave charmonium (1:4:7) and 
D-wave (1 : 1:4) []^ |^, |^. Within a charmonium picture the favoured interpretation 
has been that these ratios arise as a consequence of nodes in the radial wave function of 
'0(35) which suppress DD and DD* due to the coincidence of p (4040—> DD, DD*) lying 
near the Fourier transformed node (Figure 1) p, |^. This hypothesis is quantitatively 
consistent with the model description of other cc dynamics P and we shall therefore adopt 
this as our point of departure. We note a useful kinematical coincidence. 


p(4160 ^ DD*) 0.75 

p(4040 ^ DD) ~ (h77 ^ ^ 

p(4160 ^ _ 0.54 

p(4040 ^ DD*) ~ (h^ ^ ^ 

If the decay amplitudes are insensitive to these small changes in momenta then one would 
expect immediate correlations along the following lines. If the DD and DD* channels 
of 0(4040) are suppressed due to the node structure of the (35) wave function, and if 
the hadronic decays of both resonances 0(4040) and 0(4160) are due to a common 0(35') 
component, then one infers that there will be a corresponding suppression of DD* and 
D*D* for the 0 (4160). After taking due account of the spin weightings, one expects 


F(4160 DD*) ~ 4F(4040 DD) 
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Table 1: Width ratios in the flux-tube model [pi] ] in regions A, B and C of parameter space. 
The various width ratios are calculated at each “dot” in Figure ^ and the mean R and standard 
deviation a are computed for each region. The “dots” in Figure]^ fit the experimental ratios in eqn. 
^ up to IcT deviations given in ref. |19|. This translates into cr/i? ~ 40%. Ratios indicated with 
a/R significantly smaller than this should hence be understood to contain parameter independent 
information, and are hence more model independent. No assumption about total widths are made. 
The first ratio would still be valid for ■(/’(dOdO) as pure 3S. The decays V'- —> DgDg in Region B 
and '0+ —> D*Ds in Region C are at a node and hence very sensitive to parameters, indicated by 
an asterisk. The conventions and parameters are those of ref. |^. 
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Region A 


R 


4.0 


36 24% 
1.4 
3.3 
4.1 


2.3 

2.4 
3.9 


1.3 

15 


a/R 


20 % 


5.8 
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0/0 

3% 

3% 
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Region B 


R 


13% 

8.8 

.25 

10 


5.6 

1.2 

* 


.69 

5.0 


a/R 


29 

45% 

58% 

30% 


O/o 

5% 

* 


26% 

9% 


Region C 


R 


54 


29% 

* 

3.4 

2.9 


2.1 

1.8 

* 


1.1 

16 


a/R 


42% 


* 

48% 

8 % 


14% 

5% 

* 


23% 

24% 


r(4160 ^ D*D*) 


~ -r(4040 ^ DD* 
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from which one can dednce the naive expectation r(4160 —> D*D) : r(4160 —> D*D*) = 
1 : 3.5li;| using the ratios in ref. 0. Detailed calculations following the parameters, 
conventions and methods of ref. p[ show modifications to these numbers indicated in Table 


The positions of the nodes in momentum space depend on the spatial spread of the 
wave functions which are summarised by a parameter [3 - see Figure ||. Depending on the 
magnitude of (3 one can find different bands of solutions in (3 space that fit the branching 
ratios of V’(4040) ^ DD ; DD* ; D*D* (see Figure H). Historically the solutions^ of refs. 
|T8| , p4| appear to be those of band “A” where the outgoing momenta pd*d and pdd are 

both in the vicinity of the first node, thereby suppressing the DD and DD* channels. A 
similar situation ensues for band “C” where both momenta are in the vicinity of the second 
node. Solution A gives 


a(e+e“ DD : DD* : D*D* : D,D, : D,D*) ~ 3 : 1 : 14 : 0.4 : 10 

^ Pa = P = 0.44 GeV [j^ was fitted from experiment. Pa = 0.485 GeV, P = 0.39 GeV and Pa = 0.46 
GeV, P = 0.42 GeV were deduced from numerical wave functions consistent with spectroscopy. We 
obtain with a coulomb + hnear potential that Pa = 0.41—0.48 GeV for a string tension b = 0.18±0.02 GeV^, 
a charm quark mass nic = 1.5 ± 0.3 GeV and as = 0.4 ± 0.1. Restricting to & = 0.18 GeV^, as = 0.4 and 
rric = 1.5 — 1.8 GeV preferred by potential models §, we predict Pa = 0.45 — 0.48 GeV. {Pa, P) is defined 
in Figure 
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Solution C gives qualitatively similar results to solution A (except for the DgD* mode) 

cr(e+e" ^ DD : DD* : D*D* : D^Ds : DsD^) ~ 3 : 1 : 10 : 0.3 : 0 

In solution B, poD is near the second node and hence dramatically suppressed but pdd* 
is between the hrst and second nodes with a similar magnitude, but opposite sign, to that 
in solution A. When the magnitude of p is varied by the amounts in eqs. 0-1 the effect on 
solution B is to move DD* into the region of a local maximum and to shift the DD away 
from the node (and hence to increase the amplitude) such that the effects of the radial 
nodes are washed out; as a result the rates DD ; DD* : D*D* tend to restore towards their 
“naive” spin counting value of 1:4:7. The actual ratios depend upon the parameters chosen 
within the band (see Table |^) the most probable results being for solution B 

a(e+e- ^ DD : DD* : D*D* : D^D^ : D^D*) ~ 1 : 4 : 10 : 2 : 1 

We shall consider the consequences of each of these “mathematical” solutions, even 
though the parameters of solution A are physically preferred. 

The production of Dg and D* from the '^(4040) and '^(4160) provide further consistency 
checks. The solution A which gave a dramatic set of branching ratios for the D system 
implies an effective “spin counting” result for the Dg states arising from the common 3S 
component of the ^’(4040) and ^’(4160). This is because 

^(^^(4160) ^ DgD*) = 0.41 GeV 

^(■^(4040) — DgDg) = 0.45 GeV 

and hence their momenta are far from the nodes of the (OS') component wave function. On 
the other hand, 

^(■^(4160) ^ DgDg) = 0.67 GeV 

is coincident with the node: so in this scenario one has the result that, due to the common 
3S component, the spin counting is realised best between the two states 

r(?/>(4160) ^ DgD*) ~ 4r(V’(4040) ^ DgDg) 

while 

T{ij{m0) DgDg) ~ 0 . 

This characteristic set of results is very different from those in solution B where prox¬ 
imity of the node suppresses r('0(4O4O) —> DgDg) ~ 0. 

Thus in summary the three patterns that arise in this '0(35') — He mixing picture of the 
■0(4040) and 0(4160) are : 

Solution A (physically realistic) 

0(4040) ^ DD< DgDg < D*D < D*D* 

0(4160) ^ DgDg < D*D <DD< DgD^ ^ D*D* (6) 

Solution B 

0(4040) ^ DgDg <DD < D*D < D*D* 

0(4160) ^ DD^ DlDg ^ DgDg < D*D < D*D* (7) 
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Solution C 


^/>(4040) ^ DsD, < DD < D*D < D*D* 


^/>(4160) 


D;Ds < DsDs < D*D <DD < D*D^ 


( 8 ) 


with quantitative measures given in table 1 . 

Independent evidence for the ^’(4040) and '^(4160) consisting of '0(3*5') mixed with an 
inert state arises within this model from a study of their total widths. Figure 2 shows 
that the regions in parameter space consistent with the experimental widths are the same 
independent of the pair creation amplitude, implying a common dynamics. 

If the tests in eqs. ^ are successful they will conhrm the common presence of the (SF) 
component in the wave functions and that the additional component is inert to hadronic 
decays thereby implying a new dynamics beyond that in conventional charmonium. Proving 
that this is due to hybrid states would then require further signals in other partial waves 
since the above mixing will ironically have made the hybrid dynamics effectively invisible 
in the 1 channel. For example, vector hybrids would not then be a signihcant source of 
enhanced '0(3685) at CDF [^. However, if the exotic 1 ^ state lies below the '0(4040) 
or '0(4160), there is the possibility of a direct spin flip Ml radiative transition of the 
hybrid component He —> 7 I ^ exposing the exotic hybrid charmonium unambiguously. 
The matrix element for 7 I is related to that of '0 —>■ rje'y in the limit where the 

photon momentum tends to zero. The Ml transition involves cc spin flip from S' = 1 to 
S' = 0, the only additional feature involving the coupling of spin and the orbital angular 
momentum associated with the excited flux-tube for the He = 1 state. To the extent 
that the flux-tube excitation is transversely polarised |l|, ^ it is effectively = ± 1 . 
Hence, 


M = ('0±|/^|-ffc(l "'")) = cos9 (11,10|11) p 

where fi is the strength of the Ml transition 0 —^ r/cb and cos 6 ' ~ ^ is the (0±|ifc) mixing 
angle. Thus the relative rates scale as 

implying a branching ratio of 0(10“^). The subsequent transition 

0± ^ He{l~"^) + 7 ^ 0(3095)77 

may provide the 0(3095) as a tag, though a dedicated search at a high intensity Tan Charm 
Factory may still be required to isolate this signal. 

There may be analogous signals in the T spectrum. There is a mass shift of the T (4S') 
and the candidate T(5S')(10580) that appears qualitatively similar to those of 0(4040) and 
0(4160) which can conceivably be explained by an increase of 80 MeV in the T(5S') — T(4S') 
mass separation due to coupled channel effects 0. Futhermore, there is a small discrepancy 
between leptonic widths and theory but not at any signihcant level. Any mixing with an 
“inert” hybrid here would require further states p to be seen, albeit at a low level, since the 
highest mass state T(11020) is the sixth in the tower and is consistent with being OS', hence 
militating against strong mixings with additional degrees of freedom. Given that there will 
soon be extensive studies in the 4S peak at B-factories it may be interesting to see if there 
are any anomalous radiative decays to I”’*' Hb state arising from a possible Hb component 
mixed into the 45' peak. However, we consider the charmonium states 0(4040) — 0(4160) 
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currently to be the most likely examples of hybrid mixing in the heavy quark sector and 
we urge a dedicated study of the ■^(4160) in particular to test this hypothesis. 

We wish to thank T. Barnes, J.M. Richard and W. Toki for comments. 
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Region A 



Region B 



Region C 



Figure 1: The SF —> IF + IF amplitude A = iy f (cc ^ cu + cu) in as a function of the 

CM momentum p in GeV. The physical T = SF T MbMq/Ma — SF T, where the spin-averaged 
masses M are defined in ref. |^, and the spin-flavour factor SF is 1 (DD), 4 {D*D), 7 {D*D*), 
I {DgDs) or 2 {D*Ds). We display A at typical {Pa,/3) (in GeV) : (0.475,0.455) (Region A), 
(0.290,0.265) (Region B) and (0.245,0.215) (Region C). The amplitude is slightly modified for 
r(cc ^ cs + cs) due to a s-quark mass differing from the u-quark 
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Figure 2: Total hadronic widths of ip- = ■0(4040) and 0+ = 0(4160) fitting experiment. The 
contours enclose narrow regions with experimentally acceptable [l^] total widths, indicating la 
variations from the mean experimental width. The regions where the 0_ decay ratios fit experi¬ 
ment are indicated by the “dots”. We denote them by Regions A, B and C. 
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Total hadronic widths of 0_ and 0+ fitting experiment 

























